Astronomy & Astrophysics manuscript no. aa3009-05 


February 5, 2008 


(DOI: will be inserted by hand later) 





The orbit of the close spectroscopic binary e Lup and the 
intrinsic variability of its early B-type components 
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Abstract. We subjected 106 new high-resolution spectra of the double-lined spectroscopic close binary eLup, obtained in a 
time-span of 17 days from two different observatories, to a detailed study of orbital and intrinsic variations. We derived accurate 
values of the orbital parameters. We refined the sidereal orbital period to 4'?55970 days and the eccentricity to e = 0.277. By 
adding old radial velocities, we discovered the presence of apsidal motion with a period of the rotation of apses of about 430 
years. Such a value agrees with theoretical expectations. Additional data is needed to confirm and refine this value. Our dataset 
did not allow us to derive the orbit of the third body, which is known to orbit the close system in ~ 64 years. We present the 
secondary of e Lup as a new (3 Cephei variable, while the primary is a /3 Cephei suspect. A first detailed analysis of line-profile 
variations of both primary and secondary led to detection of one pulsation frequency near 10.36 c d~^ in the variability of the 
secondary, while no clear periodicity was found in the primary, although low-amplitude periodicities are still suspected. The 
limited accuracy and extent of our dataset did not allow any further analysis, such as mode-identification. 
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1. Introduction 

Line-profile variable early B-type stars which belong to a close 
binary system are very interesting targets for studying the ef- 
fect of tidal interactions on pulsation-mode selection and/or on 
the enhancement of pulsation-mode amplitudes. As the num- 
ber of cases observed with confirmed tidally induced modes 
is sparse (HD 177 863, De Cat & Aerts 2002; HD 209 295, 
Handler et al. 2002), two complementary surveys have been set 
up to systematically study the behaviour of non-radial pulsa- 
tions (NRPs) in early B-type stars with a close companion. The 
first survey is called SEFONO (SEarch for FOrced Non-radial 
Oscillations, Harmanec et al. 1997), which searches for line- 
profile variations (LPVs) in the early-type primaries of known, 
short-period binary systems and studies orbital forcing as a 
possible source of NRPs. The second survey has a complemen- 
tary approach (Aerts et al. 1998) and starts from a sample of 
selected main-sequence early-type B stars without Balmer line 
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emission, which are known to exhibit NRPs and which turn out 
to be part of a close binary. As high-resolution spectrographs 
allow detailed study of low-amplitude features moving through 
the line profiles, both surveys are spectroscopic surveys and 
focus on the presence of LPVs as signatures of the forced os- 
cillation modes. Aerts & Harmanec (2004) have compiled a 
catalogue of line-profile variables in close binaries. 

In the framework of such a systematic study of line-profile 
variable early B-type stars, we selected known close binary 
systems from the list of candidate /? Cephei variables result- 
ing from the high-resolution spectroscopic survey for LPV in 
09.5-B2.5 II-V stars performed by Schrijvers et al. (2002) and 
Telling et al. (2003). One of the interesting results of this sur- 
vey is detection of LPVs in 16 of the 27 known short-period 
binaries (Porb < 10'^) in the sample (Telling et al. in prepa- 
ration). In this paper we present results from the study of the 
eccentric binary e Lup. 

2. Present knowledge about e Lup 

The target eLup (HD 136504, HR5708, aaooo 
15''22'°40!87, ^2000 = -44° 41' 22'.'64, my = 3.367) 
is a double-lined spectroscopic binary (SB2) in a close 
orbit (Curtis 1909; Moore 1910; Campbell & Moore 1928; 
Buscombe & Moiiis 1960; Buscombe 8c Kennedy 1962; 
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Fig. 1. Spectra of e Lup centered at the Si III triplet near 4568 A 
obtained during nights in September 1995 and April 1996 with 
the ESO CAT/CES combination (Schrijvers et al. 2002). The 
spectra are offset for clarity. 



Thackeray 1970) and the brighter component of a visual binary 
with a period of several decades (Thackeray 1970, hereafter 
T70). Buscombe & Kennedy (1962) published the first orbital 
solution, resulting in a nearly circular orbit (e < 0.03) with an 
orbital period of approximately 0'?901 days. Regrettably, their 
radial velocities (RVs) were never published. T70 determined 
an orbital period of 4'?5597 days and an eccentricity of 
e = 0.26. No signs of cyclic changes due to ellipticity were 
detected, and no eclipses were observed. T70 reported the 
presence of the third component (separation 26". 5) with a 
possible orbital motion of 64 years. T70 provided spectral 
types for all three components, B3IV, B3V, and most likely 
A5V. Other estimates found in the literature provide only one 
type for the whole multiple system, which ranges from B3IV 
to B2IV-V (de Vaucouleurs 1957; Levato & Malaroda 1970; 
Hiltner et al. 1969; BSC, Hoffleit 1991). No other studies of 
the orbital parameters have been carried out so far. 

The HIPPARCOS parallax resulted in a distance of 155 ± 
15pc (Perryman et al. 1997), which led De Zeeuw et al. (1999) 
to conclude that e Lup is most probably a member of the UCL 
OB association, located at a distance of 140 pc. In this case, 
the estimated age of the multiple system of eLup is about 14- 
15 Myr. 

Tokovinin (1997) gives the following values of the visual 
magnitudes of the separate components of the system AaB: 
mvAa — 3.62, — 5.10. In his catalogue we also find 

estimates for the masses, based on the mass function of the 
orbital solution presented by T70: Mb = 7.64Mq, MAa = 
24.70Mq, whereby Ma = 13.24Mq and Ma = 11.46M0. 
Based on Geneva photometry, De Cat (2002) lists M = 7.1 ± 
O.6M0 and R = 5.0 ± I.ORq, assuming a single star. The 
following estimates of the effective temperature, log g, and the 
luminosity are available in the literature: Toff = 21 230 K (vis- 
ible spectrophotometry, Morossi & Malagnini 1985); Tcff = 
23 410 ± 2 510 K (spectroscopy, Sokolov 1995); T^s = 19 240 
K, logg = 3.85 (ubvy photometry, Castelli 1991); Teff = 
19275±1800K, logg = 3.89±0.17, logL = 3.5O±O.18L0 
(Geneva photometry, De Cat 2002). 



Table 1. Journal of the observations of s Lup obtained in May 
2003. The top part is the logbook of the 56 FEROS (1 1-16 May 
2003), while the bottom part gives information on the 47 obser- 
vations made with GIRAFFE (20-26 May 2003). The columns 
list the Julian date of the observations, the number of spectra 
obtained, the mean S/N-ratio, the mean integration time ex- 
pressed in seconds, and the number of spectra for which the 
profiles of primary and secondary are separated. 



Date (HJD) 


N 


S/N 


At (s) 


Nscparatc 


2452770 


18 


450 


330 





2452771 


13 


550 


390 


13 


2452772 


6 


460 


490 





2452773 


10 


570 


405 


10 


2452774 


1 


250 


600 





2452775 


8 


400 


420 





2452780 


6 


190 


490 


3 


2452781 


13 


210 


370 


10 


2452783 


16 


200 


450 


11 


2452784 


12 


185 


365 






In the literature we find different values of the projected 
rotational velocity wsini: 166 km s^^ (Bernacca & Perinotto 
1970); 142 km s'^ (Uesugi & Fukuda 1970); 170 km s"! 
(Levato & Malaroda 1970); 40 km s^^ (Slettebak et al. 1975); 
133 km s~^ (Hoffleit 1991). The higher values are obtained by 
treating the system as a single star. 

The intrinsic variability of e Lup has not yet been studied in 
detail. Schrijvers et al. (2002) detected bumps in the primary, 
possibly caused by NRP of high-degree (£ > 2). Their spec- 
tra, taken with the CAT/CES combination at ESO, La Silla, 
during nights in September 1995 and April 1996, are given 
in Fig. n and are included in the analysis presented in this 
paper Signatures of both duplicity and line-profile variability 
of the primary are visible. Both components of the spectro- 
scopic binary are located in the theoretical f3 Cephei instability 
strip (Pamyatnykh 1999) making them candidates for exhibit- 
ing /3 Cephei type pulsations. 

3. Data and data reduction 

We obtained two quasi-consecutive time-series of high- 
resolution echelle spectra for sLup, obtained from different 
observatories, in order to cover the orbit as well as possible. 

The first part of the time-series of e Lup was gathered with 
the FEROS echelle spectrograph attached to the 2.2m ESO tele- 
scope at ESO, La Silla, Chile during the nights of 1 1-16 May 
2003. The night of 15 May was mainly lost due to bad weather. 
The spectra have a resolution R ^ 48000 and cover the range 
3500 to 9200 A, divided into 39 orders. We obtained 56 spec- 
tra with an average S/N-ratio of 450 in the spectral range of 
the Si III triplet near 4568 A. The integration times were be- 
tween 5 and 13 minutes. A logbook of the 6 nights is given 
in the top part of Table ^ The data reduction was performed 
using the on-line FEROS reduction pipeline, which makes use 
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JNr. 


Observatory 


Telescope/Instrument 


Time interval 


P KVS 


KVS 


Dispersion 








(HJD- 2400000) 


compl 


comp2 


(A/mm^^) 


1 


Lick Southern Station 


0.929m/l-& 2-prism spectrograph 


17691.76-19221.59 


10 


7 


10.3-20.3 


2 


Mount Stromlo 


0.762m/3-prism spectrograph 


35297.96-36798.86 


3 


3 


36 


3 


Radcliffe 


1.9m/2-prism spectrograph 


38456.59-40369.36 


32 


18 


6.8-15.5 


4 


ESO 


CAT/CES 


49965.60-50196.60 


3 


3 


R=65000 


5 


ESO 


2.2m/FEROS 


52770.51-52775.86 


56 


25 


R=48000 


6 


SAAO 


1.9m/GIRAFFE 


52780.28-52784.61 


47 


39 


R=32000 



1: Campbell & Moore (1928); 2: Buscombe & Moniii (1960); 3: T70; 4: Schrijvers et al. (2002) 



spectrograph at the 1 .9m telescope at SAAO, Sutherland, South 
Africa. From the 7 nights of 20-26 May 2003, 3 nights were 
lost due to bad weather Integration times were between 5 and 
10 minutes. The total wavelength range covers 4400-6680 A 
and is spread over 45 orders. Due to the lower resolution the 
GIRAFFE data are of lower quality than the FEROS data. The 
average S/N-ratio we obtained was 200. A logbook of the spec- 
tra can be found in the lower part of Table ^ The GIRAFFE 
spectra were reduced by using the GIRAFFE pipeline reduc- 
tion program XSPEC2. At the top of Fig. Ha sample of reduced 
and normalised GIRAFFE spectra centered on the wavelength 
region of the Si III triplet near 4568 A are shown. 

4. The orbital motion 

T70 studied the orbital motion of the triple system e Lup from 
spectra obtained with the Radcliffe telescope during 1964- 
1970. The orbital elements he derived are listed in the first 
column of Tabled Our dataset, consisting of 103 spectra ob- 
tained from two different observatories, samples the short or- 
bital period approximately 3 times, whereby 80% of all phases 
is covered. Our short dataset, in combination with the older 
data, does not allow investigation of the triple system, due to 
the extremely poor phase coverage. 



4. 1. Radial velocities 

Fig. 2. A randomly chosen set of reduced spectra of e Lup cen- We derived the RV values of each of the two components in two 
tered at the Si III triplet near 4568 A obtained during respec- different ways: from the first normalised velocity moment (v) 
tively 6 and 4 nights in May 2003 with the FEROS (bottom) (e.g. Aerts et al. 1992) and by means of the cross-correlation 
and GIRAFFE (top) spectrograph. The spectra are offset for (CC) technique. 

clarity. To calculate of the {v) of the three lines of the Si III triplet at 

4552.654 A, 4567.872 A and 4574.777 A, we used variable in- 
tegration boundaries. Values of the RVs of the secondary could 
only be derived at phases near elongation, when both profiles 

of the ESO-MIDAS software package. We performed an ad- were separated. One has to keep in mind that the estimates of 

ditional correction for the wavelength sensitivity of the shape the RVs of the primary near conjunction are contaminated by 

of the internal flatfields by means of a smoothed average of the presence of the secondary. 

dome flatfields. For each of the 6 nights, a randomly chosen set We appHed the CC technique to several wavelength re- 
of normalised FEROS spectra centered at the Si III triplet near gions with well-defined absorption lines for primary and sec- 
4568 A is given at the bottom of Fig. |2l ondary. These regions were centered on the following 6 lines: 

The second part of the time-series of eLup consists of 47 Sim 4568 A, Mgll 4481 A, He I 5016 A, He I 4917 A, 
spectra (R ~ 32 000) measured with the GIRAFFE echelle He i 5876 A, and He i 6678 A. As templates for the CC, we 



HJD 2452784 




4550 4560 4570 

Wavelength (A) 
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Table 3. Orbital elements and their standard errors of the close binary system of eLup, derived with the FOTELcode from the 
RVs of the individual dataset of Lick (second column) and Radcliffe (third column) data, of the combined Lick & Radcliffe data 
(fourth column), and of FEROS & GIRAFFE data (last column). Weights according to the formula w — i?/32000 were assigned 
and data points with errors larger than 10 km s^^ were removed. All epochs are given in HJD-2400000. As comparison, the first 
column gives the solution of the system, as presented by T70, which is based on 18 Radcliffe spectra. The dataset numbers are 
as indicated in Table |2| 



Nr. 3 1 3 1+3 5+6 

Elem. T70 FOTEL FOTEL FOTEL FOTEL 



Porb (days) 


4.559783 ± 2 10"" 


4.5597 fixed 


4.5597 ± 1 10"* 


4.55977 ± 1 lO"'' 


4.560 ± 8 10"^ 


To 


39370.68 ± 0.09 


18304.4 ± 0.3 


39370.71 ± 0.08 


39370.68 ± 0.08 


52767.58 ± 0.01 


e 


0.26 ±0.03 


0.26 fixed 


0.26 ±0.03 


0.26 ±0.03 


0.300 ± 0.006 


Ki (km s"^) 


56.1 ±1.5 


55 ±5 


51 ± 2 


52 ±2 


54.0 ± 0.3 


K2 (km s-^) 


64.8 ± 1.8 


71 ±5 


69 ±2 


69 ±2 


64.0 ± 0.3 




330 ± 10 


319 ±25 


332 ±8 


329 ±7 


24 ± 1 


ai sin i (A.U.) 


0.0227 


0.022 ± 2 10"^ 


0.0206 ± 8 lO""* 


0.0210 ± 8 10"* 


0.0215 ± 5 10"'' 


a2 sin i (A.U.) 


0.0262 


0.0289 ± 2 10"^ 


0.0279 ± 8 lO""* 


0.0278 ± 8 10^* 


0.026 ± 5 10"^ 


/(M2) (M©) 


0.06 ±0.01 


0.07 ±0.02 


0.056 ± 0.007 


0.060 ± 0.007 


0.065 ± 0.005 


/(Ml) (Mq) 


0.14 ±0.02 


0.16 ±0.03 


0.14 ±0.01 


0.14 ±0.01 


0.107 ±0.007 


v-y- (km s"""^) 


8±1 


j=I:5±3 


j=3: 9 ± 1 


J=l: 6 ± 3 


j=5: 3.2 ± 0.3 








j=3: 8 ± 1 


j=6: 0.4 ± 0.3 


rmsj (km s~^) 




j=I:8.8 


J=3: 7.6 


j=l: 9.0 


j=5: 2.6 










j=3: 7.6 


j=6: 2.1 


rms (km s~^) 


3.6 


8.8 


7.6 


8.1 


2.4 


-Porb.widc (yrs) 


64 











took the average normalised spectrum of the FEROS and of the 
GIRAFFE datasets. The position of the line center, measured 
by means of a Gaussian fit, was used as an estimate of the RV. 
The CC technique also allowed us to determine the RVs of the 
secondary when its profiles were semi-detached from those of 
the primary. Finally, we calculated the median RV of the 6 re- 
gions and its standard deviation. 

It turned out that the scatter on the first velocity moment 
calculated from the Si III profiles of the FEROS spectra was 
less than for the RVs calculated by CC. This is possibly due 
to the good quality and high-resolution of the FEROS spectra 
and the nicely defined absorption profiles of the Si III triplet. 
Therefore we used the values of the average of the first mo- 
ments of the 4553 A and 4568 A profiles as RVs of the FEROS 
time-series in the subsequent analysis. The GIRAFFE spectra, 
on the other hand, are of lower quality than the FEROS spectra. 
For these spectra the noise was suppressed by the CC proce- 
dure, so that we used the median of the CC measurements as 
RVs of the GIRAFFE time-series. 

Additionally, we included the three high-quality 
(R~65 000, S/N > 450) CAT spectra, obtained by Schrijvers 
et al. (2002) in 1995-1996. We calculated the first normahsed 
velocity moment {v) for both primary and secondary. 

We enlarged our sample of RV measurements with the 
ones T70 used for his analysis. These datapoints include 10 
values calculated from spectrographic plates obtained from 
Lick Southern Station, Chile, between 1907 and 1911, and 
32 RVs calculated from several spectral lines measured with 
the Radcliffe 74-inch reflector at Radcliffe Observatory dur- 
ing 1964-1970. We recalculated more accurate HJD values for 



the Lick data, based on the information given by Campbell & 
Moore (1928) and discovered that in their ambiguous notation 
the RVs of primary and secondary of nights HJD 2417703 and 
HJD 2418760 were exchanged. For 25 of the 42 datapoints, 
RVs for both components are given. Also three RV measure- 
ments (1955-1956), provided with errors, are available from 
Buscombe & Morris (1960). A journal of all available RVs we 
used for the further analysis is given in Table|2] while the indi- 
vidual RVs with corresponding HJDs can be found in Table|6| . 
The error estimates, given in Table|6lonly for our observations, 
are the rms errors of the mean of Si III, Mg II, and He I RVs. 

4.2. Determination ofttie orbital parameters 

Analysis of RVs was carried out using the FOTELcode 
(Hadrava 1990), which is designed to solve the light- and/or 
RV curves of binary stars with a possible third component and 
also can model secular changes in the orbital period and apsidal 
motion. To each individual RV we assigned a weight which is 
proportional to the spectral resolution R: w — i?/32000. The 
spectral resolution can be expressed by the formula 



with A the central wavelength of the spectrum in question in 
A, D its linear dispersion in A/mm^^, s the "pixel" spacing in 
mm, and n the number of dA per FWHM of the projected slit 
width, which has a typical value of 2-3. We adopted the values 
5=0.020 mm and n=2 for the photographic spectra. Information 

' Table|6|is only available in electronic form. 
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Nr 


1 ^ f\ 


3-6 








FOTEL 


FOTEL 




KOREL 


Jorb.ano l,1dys; 


A ^^QR'^ -1-11 


4.5598 ± 3 10" 


-4 


A ^^QR'X fiYpH 
4:.uOyoO IIACU 


Jorb.sid (.aays; 


A '^f^Qvn -1-11 


4.5597 ± 3 10" 


-4 







QQQ7n 7R -1- n 

oyo 1 u. 1 o zc u.uo 


52790.35 ± 0.02 


oz / yu.oo 


6 


n 979 _|_ n nn/:: 


0.271 ± 0.005 




977 

U.Z ( ( 


XV 1 l^KIIl S ) 




53.9 ±0.3 




Oo.o 






64.0 ± 0.3 






(j; ( ) 


0/17 _|_ c 


18 ±2 




1 7 


ui { lyi) 


n Q _i_ n 9 
u.o zn u.Zi 


1.1 ±0.2 




u.o 11 ACQ 




400 _|_ 7Q 

4:ZrO zn ( y 


319 ± 55 




4:ZO llACU. 


(Al Dill t • / 


n noi 7 + 91 n-4 


0.0217 ± 1 10" 


-4 


0.0216 


n n QAn /J lA 111 


n n9c;e -X- 9 1 n~'* 

U-UZiJO ZL _LU 


0.0258 ± 1 10" 


-4 


u.u^u± 


J ) \A^J^Q) 




0.066 ± 0.001 




0.065 


J{lVll) \1V1Q) 


n 1 1 4- n 009 


0.110 ±0.002 




114- 

U. -L J-^ 


1 1 /'t'Tn c ^\ 
V'y-^ V^Klll a ) 


^ m 








1 1 /^Irm c ^\ 


9 ± 2 


y zt z 






1 1 /^ttn c ^\ 


m 


3zt3 








1 1 -L n 9 

J L m u.^ 


1.1 zb 0.2 






1 1 f^Irm c ^ 


— + 
u.o zn u.o 


-3.0 zb 0.2 






nn c -1 f^Irm c A 
llllai l^R-lil a ) 


11.2 








rmss (km s^^) 


9.8 


8.2 







rms4 (km s~^) 


5.8 


5.3 






rmss (km s~^) 


1.9 


1.8 






rmse (km s~^) 


2.2 


2.2 






rms (km s~^) 


2.9 


2.6 







Table 4. Orbital elements and their stan- 
dard errors for the close binary system of 
eLup, derived with the FOTEL code from 
all RV data, except Mt Stromlo (first col- 
umn), and from the subset of Radcliffe, CAT, 
FEROS, and GIRAFFE data (second col- 
umn), including apsidal advance. In the pre- 
sented solutions, the RVs of CAT, FEROS, and 
GIRAFFE were replaced by the correspond- 
ing KOREL RVs. Weights according to the for- 
mula w = i?/32000 were assigned and data 
points with errors larger than 10 km s^^ were 
removed. The last column gives the results as- 
sociated to the best KOREL disentangling solu- 
tion of CAT, FEROS, and GIRAFFE spectra, 
in which the values of Poib and lo were kept 
fixed. Both anomalistic and sidereal period are 
hsted. All epochs are given in HJD-2400000. 
The dataset numbers correspond to the ones in 
TableEl 



on the dispersion or resolution of the different types of data is 
provided in Table |2] 

First, we calculated a local solution for the subsets of 
Radcliffe, Lick, and the combined FEROS and GIRAFFE data. 
In the case of the Lick data we kept the values of the orbital 
period and eccentricity fixed. For the new data, measurements 
with errors larger than 10 km s^^ were rejected. These uncer- 
tain measurements correspond to datapoints near conjunctions. 
The lack of accurate datapoints close to the system velocity 
due to blending of the lines is a typical problem that occurs 
for SB2s. One of the consequences is that orbital parameters 
whose determination is sensitive to the velocity behaviour near 
conjunction, such as the eccentricity e and the periastron longi- 
tude uj, will be less accurately derived. The resulting param- 
eters of the sub-datasets are listed in Table |3l From the old 
dataset of Lick and Radcliffe measurements we found a so- 
lution consistent with the one T70 derived from the Radcliffe 
data (Porb = 4'?5597, e = 0.26). The RVs of the FEROS and 
GIRAFFE time-series could be folded with the same orbital 
period with a slightly higher eccentricity e = 0.311. We note 
that the values of lo differ in the different datasets and hence 
suspected the presence of an apsidal advance. 

Secondly, we searched for solutions in larger subsets, 
thereby allowing for apsidal advance. We worked iteratively 
in cases of subsets including datasets 4-6: we determined val- 
ues of Porb and cj, which were subsequently kept fixed in the 
KOREL disentangling process of all FEROS, GIRAFFE, and 
CAT spectra (see Sect. 14. 3t . As KOREL allows construction of 
a precise RV curve for both components, even at phases near 
conjunction, we subsequently ran the FOTEL code again, re- 



placing the FEROS, GIRAFFE, and CAT measurements by the 
KOREL RVs to optimise the solution. 

We investigated the rate of the rotation of the line of ap- 
sides in the subsets including Lick and Radcliffe data (datasets 
1 and 3, third column TableQ and Radcliffe, CAT, FEROS, and 
GIRAFFE data (datasets 3-6, using KOREL RVs for datasets 
4-6, second column Table |4}. The bad phase coverage of the 
Lick data, in combination with its poor quality, did not allow 
good definition of lo, so we did not find a significant apsidal 
advance in the first subset. Next, we combined all datasets, 
spanning some 96 years, including the few Mount Stromlo 
and CAT velocities, and searched for the best matching or- 
bital solution. We did not succeed in including the few RV 
measurements published by Buscombe & Morris (1960), as 
they fell completely outside our velocity curve. Therefore we 
question the correctness of these values. The more so because 
Buscombe & Kennedy (1962) published an unreliable orbital 
solution with a nearly circular orbit of 0.9 days. Leaving out 
the Mount Stromlo measurements led to an orbital period of 
P^^^ = 4':i55983 and w = 0.8 ± 0.2°/yr, the latter value 
being slightly smaller than obtained from the datasets 3-6. 
Including the corresponding KOREL RVs for CAT, FEROS, and 
GIRAFFE data led to the 'best' orbital parameters, i.e. with 
the lowest rms value, given in the first column of Table |3 
We stress, however, that although the solutions are quite sta- 
ble concerning parameters Porb, e, Ki, and K2, uncertainties 
surround the values of lo, lo and consequently also Tq. In order 
to gain clarity, additional datasets, strategically chosen in time, 
are needed. 
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The phase diagrams of Lick & RadcHffe data (4th column, 
Table|3}, FEROS, & GIRAFFE data (last column, TableQ and 
FEROS, GIRAFFE, and CAT data (KOR EL solution, last col- 
umn, Table|4} with respect to Poib displayed in Fig.|3] 

The system velocities of the newly obtained data differ sig- 
nificantly from the older dataset 3 (see Table0}, which supports 
the presence of a third body. Given the solution of the close bi- 
nary system as derived from the whole available dataset, we at- 
tempted to converge to a solution of the triple system, but - not 
surprisingly, given the poor phase coverage - without satisfac- 
tory results. To unravel the triple system with an orbital period 
of approximately 64 years, a dedicated long-term project is re- 
quired with the aim of gathering datapoints that are well-spread 
in the orbital phase. 

4.3. Spectral disentangling 

An advantage of the KORELcode (Hadrava 1995) is that it al- 
lows the study of a spectral range as a whole and thereby can 
extract information from the blended lines. We applied the KO- 
REL disentangling technique to the set of CAT, FEROS, and 
GIRAFFE spectra and focussed on the wavelength region cen- 
tered on the two bluest lines of the Si III triplet. One spectrum 
was removed from the dataset due to its poor quality. We in- 
cluded apsidal motion to account for the time gap of 2574 days 
between the CAT data and the data taken in May 2003. We 
note, however, that given the small amount of CAT spectra, 
the solutions allowing for apsidal advance do not differ signif- 
icantly from the ones not taking apsidal motion into account. 
We kept the orbital period fixed on its best FOT EL solution 
(^orb.ano = 4^55983) and searched for the most satisfactory 
solution of orbital elements in combination with the 'best' de- 
composition of the spectra. The resulting elements are listed in 
the last column of Table|3 and the derived RVs of primary and 
secondary are plotted versus orbital phase in the lower panel of 
Fig-El The normalised disentangled spectra centered at 4560 A 
are given in Fig.|3 



5. Limits on pliysicai eiements 

5. 1. Estimates of basic physical elements 

An important unknown in the unravelling of the physical pa- 
rameters of binary components is the orbital inclination, which 
can only be accurately determined in case the star is an eclips- 
ing binary. In the Uterature no evidence for the presence of 
eclipses in eLup is found. We looked at the HIPPARCOS 
lightcurve of e Lup. The data with quality label are plotted 
versus orbital phase in Fig.|5l We see that one datapoint shows 
a decrease in brightness. This happens to be the measurement 
with by far the largest error bar, which is moreover measured 1 5 
minutes apart from the brighest datapoint in the dataset, span- 
ning more than 1000 days. Therefore we consider this insuffi- 
cient evidence to claim the presence of an eclipse. Assuming 
that eLup is a non-eclipsing system implies an orbital incli- 
nation lower than 50 degrees (lorb < 50°), given the B-type 
nature of the two components. 



Table 5. Representative values of logTcff, radius, and mass 
for the components of e Lup (indicated by asterisks in Fig.|6j, 
and the associated estimate of the orbital inclination and semi- 
major axis a. 



star 


logTeff 


R/Rq 


M/Mq 


iorb (°) a (Rq) 


1 


4.34 


4.7 


8.7 


20.5 29.2 


2 


4.31 


3.8 


7.3 





We estimated the radius R of the primary by using the for- 
mula 

log(i?/i?0) = 7.474-2 logTeff-0.2BC- Ay-log TT, (2) 

with TT the HIPPARCOS parallax (6.06 ± 0.82 mas, Perryman 
et al. 1997), BC the bolometric correction, reddening Ay, 
effective temperature Tcs, and V magnitude. From the 
HIPPARCOS photometry transformed to the Johnson V mag- 
nitude by using Harmanec's (1998) transformation formula and 
from the Geneva photometry, we find V = 3™375 and V = 
3"?369, respectively. The corresponding dereddened value of 
the latter is Vq =3™32. Next, we considered two situations: that 
the two components are equally bright in the V magnitude, im- 
plying Vb,i =4™06, and that the secondary is half as bright as 
the primary, hence Vb.i =3"?75. These two situations are re- 
alistic limits on the true value of Vb,i- Considering a range of 
log Toff.! G [4.280; 4.364], which is consistent with a spectral 
type between B3 and B2 according to Harmanec's tabulation 
(1988), formulalJlleads to a radius of the primary between 3.6 
and 6.7 Rq (black box in Fig.|6jl. 

By comparing our observations with Schaller's et al. (1992) 
and Claret's (2004) evolutionary models (see Fig. |6|l, and as- 
suming an age of 15 Myr (dotted line) of the eLup system 
as a member of the UCL OB association (De Zeeuw et al. 
1999), we find the following restrictions on the radius, mass 
and log Toff of the primary of eLup: Mi e [7.6, 9.7]Mo, 
logTcff,! e [4.32; 4.36] and i?i e [4.1,5.6]i?0 (dark gray 
box in Fig. |6j. Using this range of masses for the primary 
we can narrow down the interval of the orbital inclination and 
of the mass of secondary. From the values Mi sin^ Zorb ~ 
O.37A/0 and A/2 sin"^ iorb ~ O.31A/0, derived from the 
best FOT EL solution, we find i e [19.5; 21.5]° and M2 S 
[6.4; 8.2]A/0. These values of Ali sin"^ i are thus compatible 
with the lack of observational evidence for eclipses. The evo- 
lutionary tracks (see Fig. |6jl for a secondary with M2 G 
[6.4; 8.2]A/0 and the same age as the primary, impose the fol- 
lowing limits: logTeff,2 G [4.28; 4.33] and i?2 G [3.5,4.4]i?Q 
(light gray box in Fig. |6|l. The physical parameters of repre- 
sentatives of the primary and secondary of eLup, indicated by 
asterisks in Fig.|6]and given in Table|5] are compatible with a 
B3IV and B3V star as suggested by T70. 

From the Full Width Half Intensity (FWHI) of the disentan- 
gled profiles (Fig. 13, we measured projected rotational veloci- 
ties according to the formula vi sinzj-ot — FWHI/Vlnl6, and 
arrive at values vi sinirot — 42 km s^^ and V2 sinijot = 32 
km s^^ for the primary and secondary, respectively. These val- 
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ues are compatible with the value of 40 km s^^ Slettebak et al. 
(1975) listed in their catalogue. Assuming Zj-ot G [19.5; 21.5]° 
leads to the following limits on the equatorial rotational ve- 
locities: vi € [115; 125] km s^^ and V2 e [88; 95] km s"^ 
For a primary and secondary with radii as derived above, 
we estimate a rotational period between 1^7 and 2*? 7 days 
(/rot G [0.41; 0.60] c d-i), respectively, Y}9 and 2':'5 days 
(/rot € [0.39; 0.54] c d~^). Given all uncertainties involved, it 
is quite possible that both components have identicaJ rotational 
periods, possibly corresponding to the spin-orbit synchroniza- 
tion at periastron, observed for many eccentric -orbit binaries. 
For e Lup its value is 2'?481 if the eccentricity e — 0.277 from 
the KOREL solution is adopted (cf. formula 7 and relevant dis- 
cussion in Harmanec 1988). 

We have to bear in mind that, due to several limitations on 
our dataset, the derived values of Mi, Ri, Tcs.i (i = 1,2) and 
the orbital inclination are only rough estimates, derived under 
certain assumptions. Additional information from future inter- 
ferometric observations would provide a much more accurate 
determination of the component masses and radii, as well as 
the orbital inclination and the relative luminosities of the com- 
ponents. 

5.2. Interpretation of the observed apsidal motion 

The range of apsidal periods we derived in Sect. 14.11 {U = 
428 ± 79) is consistent with the values obtained for other bi- 
naries with orbital periods between 4 and 5 days (Petrova & 
Orlov 1999). As a comparison, we mention the system QXCar, 
which is very similar to e Lup, and consists of two B2 V stars 
orbiting in an eccentric (e = 0.28) orbit of 4.5 days, which has 
an apsidal period of 370 years. 

The tide-generating potential of the binary system can give 
rise to a periodic variation of the longitude of periastron. It 
is well-known that the observed apsidal advance in eccentric- 
orbit binaries arises from both classical and general-relativistic 
terms, i.e. 

Wobs — ^cias + ^rci- We Calculated the contribution of the rela- 
tivistic term from the formula given by Gimenez (1985), which 
resulted in Wroi G [0.113; 0.139]°/yr and which is a 6-7 times 
smaller value than our observed range of values. 

An estimate of the total contribution of classical effects can 
be derived from uidas = ^obs ^ f-^rci and leads to ujcias G 
[0.687; 0.661]°/yr. The classical effects are caused by tidal and 
rotational distortions, and, if relevant, the presence of a third 
body: Wcias = t^tidai + t^rotation + ^^3- The latter term is neg- 
ligible in the case of eLup due to the long orbital period of 
the triple system. According to the formula given by Wolf et 
al. (1999), the contribution of the third component to the rota- 
tion of the line of apses is estimated as = 3 x 10"®°/yr. 
Neglecting this effect, we estimated a mean apsidal constant k2 
(cf. formula 4 in Claret & Willems 2002) for the primary and 
secondary from the observed apsidal period and relative photo- 
metric radii of the components, ri = 0.16 and r2 = 0.13, as 
k-z = 0.0074 for the representative set of values from Table|5] 
This can be compared with the theoretically expected values, 
interpolated from Claret's (2004) models, k2 = 0.0068 and 



k2 — 0.0074 for the primary and secondary, respectively. It 
is seen that there is no obvious conflict between the observed 
and theoretically predicted values. Given our observed range 
of primary radii, assuming Zorb = 20.5°, e = 0.277 and 
k2 = 0.007, we estimate the rate of secular apsidal motion 
due to the tidal distortion of the components from the for- 
mula given by Claret & Willems (2002, their formula 13): 
0.004°/yr < wtidai < 0.019°/yr As we could derive only 
rough estimates of the basic physical properties of both stars, 
we do not take into account the possible effect of dynamic tides 
studied by Claret & Willems (2002). We note, however, that 
their effect will probably be small as the stars are rotating close 
to spin-orbit synchronization at periastron unless there is some 
resonance present in the system. This must be left unanswered 
until much more accurate masses and radii of components are 
known. The estimated values of a;tidai and uj^ lead us to con- 
clude that the term Wiotation, accounting for rotational effects, 
must contribute significantly to the apsidal motion. 

6. Intrinsic variability 

High-degree (£ > 2) modes are suspected in the primary of 
eLup (Schrijvers et al. 2002; Fig. 0. The analysis of LPVs 
of both primary and secondary is a challenging task due to the 
comparable line-strength of the two components and the blend- 
ing of the component lines in phases far from elongation. The 
spectra of two line-profile variable components that have been 
disentangled by KOREL have never been the subject of an LPV 
study before. We investigate the two bluest Si III profiles be- 
fore and after KOREL disentangling and then compare the re- 
sults. A disadvantage of using the original spectra is that only 
spectra obtained at phases of elongation can be used, while 
all KOREL disentangled profiles can be considered for analy- 
sis. However we question whether the intrinsic variations of 
primary and secondary are properly disentangled, as KOREL is 
not constructed to treat pulsational variations. For the testcase 
K Sco, Harmanec et al. (2004) and Uytterhoeven et al. (2005) 
obtained good results concerning the retrieval of intrinsic vari- 
ations with high amplitude after disentangling but the perfor- 
mance of KOREL was less clear for low amplitude variations. 

The set of disentangled profiles of the first (second) com- 
ponent were constructed by adding the KOREL residuals calcu- 
lated in the restframe of the primary (secondary) to the nor- 
malised disentangled profile of the primary (secondary). The 
'original profiles' of the first (second) component were ob- 
tained from the observed spectra by a shift in wavelength ac- 
cording to its corresponding orbital velocity in the best-fitting 
orbit obtained with KOREL (last column in Table |3 bottom 
panel in Fig. O. In Fig. we compare the original and dis- 
entangled profiles of the primary and the secondary at three 
different orbital phases. From these comparison plots we con- 
clude the following: 1 . The primary and secondary of e Lup are 
line-profile variables. 2. The KOREL residuals in the restframe 
of the primary (secondary) contain signatures of the intrinsic 
variability of both components. 

In order look closely for signs of intrinsic variability, 
we show the grayscale representations of the disentangled 
Sim 4553 A profiles of the primary and secondary on two 
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Fig. 8. Grayscale representations of the disentangled FEROS spectra of the primary (top) and secondary (bottom) obtained during 
the nights of HJD 245277 1 (left) and HJD 2452773 (right). At the bottom of each figure the corresponding nightly average profile 
is plotted as a reference. 



nights during the FEROS observing run when the two compo- 
nents are near elongation (Fig.|8}. To bring out the variability 
of the primary and secondary in full detail, we calculated the 
residuals with respect to the nightly average spectrum. For both 
components we find signatures of intrinsic variability, although 
weak in case of the secondary, in the form of moving black and 
white bands. 

We performed CLEAN (Roberts et al. 1987), Lomb-Scargle 
(Scargle 1982) and Phase Dispersion Minimisation (PDM, 
Stellingwerf 1978) analyses on the first three velocity moments 
((u), (u^), (v^)), the EW, and the RVs associated to the mini- 
mum of each line-profile, in search for intrinsic variations of 
primary and secondary. Also the 2D-analysis Intensity Period 



Search (IPS, Telting & Schrijvers 1997), to find variable sig- 
nals in individual wavelength bins (AA = 0.02 A in case of 
eLup), was performed by using the CLEAN algorithm. All 
the diagnostics were calculated from the Si III 4552.654 A and 
4567.872 A profiles of both the disentangled and original pro- 
files. 

We investigated the variable signal of the complete 
datasets, including all datapoints, and of reduced datasets, 
which included only the 56 datapoints obtained when the lines 
of the component spectra were separate (for a logbook, see 
Table The clearest signals were detected in the variability 
of the latter datasets for the obvious reason that unreliable dat- 
apoints near conjunction were not included. In theory, all dat- 
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apoints after disentangling should contain reliable information 
about the variability of the primary (secondary), but the obser- 
vations show that even for the disentangled profiles, this re- 
striction to the 56 separate spectra is to be preferred. For the 
moment it is unclear how well the KOREL disentangling proce- 
dure can separate the intrinsic variability contributions of the 
primary and secondary near conjunction. This can only be de- 
duced from a dedicated simulation study, which is beyond the 
scope of the current work. 

Below, we report on the results obtained from the reduced 
datasets containing 56 datapoints and spanning 12 days. The 
window function is very complicated and leads to severe alias- 
ing. The half width at half maximum of the amplitude of the 
highest window peak is as large as 0.03 c d^^ and gives a 
rough estimate of the error on the frequencies. According to 
the empirically derived formula by Cuypers (1987), the small 
dataset allows a frequency accuracy of 0.01 c d^^. Hence, we 
searched for frequencies in the frequency domain between 
and 20 c d~^, with frequency steps of 0.01 c d^^. A more ex- 
tended discussion of the frequency analysis results than pre- 
sented here can be found in Uytterhoeven (2004)^. 

6. 1. Intrinsic variability of the primary 

No clear dominant frequency is found in the variability of the 
primary. The highest amplitude is found between 0-2 c d^^. 
In the Scargle periodogram of the non-prewhitened data, as 
shown in the upper panel of Fig.|9j one sees hints of the pres- 
ence of a variability near 6.46 c d^^ or one of its aliases, but 
these peaks reduced in amplitude (original profiles) or almost 
completely vanished (disentangled profiles) after prewhitening 
with the highest amplitude low frequency peak (lower panel of 
Fig.|9}. Signatures of the peak 6.46 c d^^ are also found at very 
low amplitude in the IPS analysis. 

The EW of the primary's original profiles shows variations 
with frequency 0.83 c d^^. This value is also recovered in the 
second moment and in the IPS analysis, and can be related to 
0.44 c d^^, detected in the EW variations of the disentangled 
spectra, through a peak of the window function. As stellar ro- 
tational frequencies are not uncommonly detected in the EW 
of line profiles (e.g. k Sco, Uytterhoeven et al. 2005; V2052 
Oph, Neiner et al. 2003; (3 Cephei, Telting et al. 1997), we 
might consider this frequency to be related to the rotation of 
the star. Indeed, the value 0.44 c d^^ is compatible with the ro- 
tational period assuming spin-orbit synchronization, as derived 
in Sect.0 

6.2. Intrinsic variability of the secondary 

The line profiles of the secondary of e Lup show variations with 
a frequency near 10.36 cd~^ (hereafter called /i) and with a 
velocity amplitude of 1.0 ± 0.2 km s~^. This frequency, or an 
alias, appears in all diagnostics after prewhitening with 0.39 
c d^^ or one of its aliases (Fig.llO>. A phase diagram of the first 
velocity moment of the disentangled profiles folded with 10.36 

^ PhD thesis available from 
ihttp://www.ster.kuleuven.be/pub/uytterhoeven_phd/| 



c d^^, is given in Fig.^] The frequency 10.36 c d^^ is not 
recovered from the IPS analysis. The EW of the secondary's 
original profiles varies with 0.75 c d^^, or one of its aliases. 

We note that the frequency 0.05 c d^^, detected in the IPS 
analysis of the profiles of the secondary, also appears in the 
IPS analysis performed on the profiles of the primary. This fre- 
quency is most probably a reflection of the orbital motion in 
the profiles, given the following relation within the frequency 
resolution: 0.05 ^ 1/4/oib- Similar relations between intrin- 
sic frequencies and the orbital frequency are observed in other 
close binary stars (e.g. a Vir, / = 1/12/orbi Smith 1985; a 
Sco, / = l/4/o,.b, Chapellier & Valtier 1992). 

The period l/fi of about 2 hours is relatively short for 
(3 Cephei type pulsations. A few examples of /3 Cephei stars 
with similar frequency values are known (for an overview see 
Aerts & De Cat 2003). So far, the highest frequency reported 
in the literature amounts to 15 c d~^ and corresponds to a very 
low RV amplitude (lu^ Sco, Telting & Schrijvers 1998). The 
low-amplitude frequency we found in the secondary of e Lup 
resembles the low-amplitude frequencies /{ = 10.48 c d~^ 
and = 10.73 c d^^ detected in the B-type primary of the 
spectroscopic binary tp'^ Orionis (Telting et al. 2001), which 
were interpreted in terms of high-degree pulsation modes. To 
confirm the presence of /i and refine its value, additional mon- 
itoring of e Lup is required. As the intrinsic variability of the 
secondary can only be studied at orbital phases near elonga- 
tion, one should aim to gather data during these phases. 

Frequency /i shows up in both original and disentangled 
datasets. This result indicates that the variability inherent to the 
secondary is well preserved through KOREL disentangling dur- 
ing orbital phases near elongation. 

7. Summary and discussion 

From a dataset of 106 high-resolution spectra obtained quasi- 
consecutively from two different observatories with a time- 
span of 20 days, we were able to refine the orbital parameters 
of the close orbit of the triple system e Lup by using the FO- 
TEL and KOREL codes. As the star is an SB2, we obtained RVs 
for both components and found an eccentric (e = 0.277) orbit 
of 4'?55970 days. By adding our data to published RV mea- 
surements we found strong evidence of the presence of apsidal 
motion (U ~ 430 years). In order to solve the triple system, 
as well as to further investigate the effect of apsidal motion, an 
extensive dataset, spanning several decades, is required. 

Two stars of spectral type between B3 and B2 in the close 
orbit, as suggested before by T70, agree with our orbital solu- 
tion. We estimated the component masses Mi e [7.6, 9.7]^^© 
and M2 e [6.4; 8.2]Mq. Both components may be spin-orbit 
synchronized. Interferometric measurements would be a mayor 
step foreward in accurately determining the physical parame- 
ters of the system. The triple system e Lup is a very suitable 
target for large interferometers. 

We used the KOREL disentangling technique as an interme- 
diate step in the study of the two early-B type line-profile vari- 
able stars of similar brightness. The variability of primary and 
secondary seems to be preserved well after disentangling. 
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Although the grayscale representations clearly show signs 
of the presence of variability, we were not able to detect a domi- 
nant frequency in the line profiles of the primary. We thus clas- 
sify the first component of e Lup as a /3 Cephei suspect. On 
the other hand, a candidate pulsational frequency was found in 
the variability of the secondary, near 10.36 c d^^. Its precise 
value could not be determined due to the poor frequency res- 
olution. This frequency certainly needs confirmation by means 
of additional intensive monitoring. Nevertheless, we propose 
the secondary of e Lup as a new (3 Cephei variable. 

Further investigation of the intrinsic variability of both 
components of e Lup requires a follow-up multi-site campaign 
with intensive monitoring of the star, preferably at phases near 
elongation, by means of high-resolution spectroscopy. For an 
example of the benefit of such a multi-site campaign, we refer 
to Handler et al. (2004) and Aerts et al. (2004b). 

The detection of several observational cases of tidally in- 
duced oscillations is necessary to improve our understand- 
ing of the relation between tidal interaction and the excitation 
of pulsation modes, and the seismic modelling of such stars. 
Effects such as the deformation of the spherical shape of a 
star due to tidal forces are usually not taken into account in 
theoretical calculations. Aerts et al. (2002) have shown that, 
in the case of the 5 Scuti star XX Pyx, deformation of the 
star due to the tide-generating potential is clearly more impor- 
tant for interpreting the pulsational behaviour than is the ro- 
tational deformation. The tide-generating potential is propor- 
tional to eT ~ {Ri/af'{M2/Mi), while the deformation due 
to centrifugal forces is proportional to ec — (/iot//puis)^- 
The values derived by Aerts et al. (2002) for XX Pyx were 
ej^ = 6 X lO^'^ and ec = 5 x 10^'*. We also calculated the 
contributions of the tidal forces and rotation on the deviation 
of spherical symmetry for the secondary of e Lup. Assuming 
iorb = 20.5°, i?i ^ 4.7i?0 (Table|5} and assuming spin-orbit 
synchronization /lot = 0.403 c d~^, we derived et ^ 0.004 
and ec — 0.002; hence, ex is larger than ec- As a compar- 
ison we calculated ec and ex for the hybrid 7Doradus star 
HD 209 295 for which tidally induced oscillations indeed have 
been detected (Handler et al. 2002). As not all parameters are 
given by these authors, we can only give a rough estimate: the 
value of Et is on the order of 3 x 10^'^, while ec is 3 orders 
of magnitude larger. For this system, the rotational effects on 
the oscillations are thus larger than the tidal effects. These ex- 
amples show that variety of behaviour exists amongst close 
binary systems. Only confrontation between several observa- 
tional findings and theory can bring us closer to understanding 
the mutual interaction between the mechanisms that play a role 
in the properties of oscillations in close binary systems. 
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Fig. 3. Observed RVs of e Lup derived from the FEROS (•) and 
GIRAFFE (o) data (top) and older Lick (biack square) and Radcliffe 
(open square) data (middle) versus orbital phase (Porb = 4^560 and 
P^^^ = 4^55977, respectively). The dark gray full and dashed 
lines represent the local orbital solutions obtained with FO- 
TEL (last two of columns Table[3j for the primary (black) and 
secondary (gray), respectively. Note the different shapes of the 
RV curves, reflecting the slow apsidal advance. In the bottom 
figure the RVs calculated from the CAT (★), FEROS (•), and 
GIRAFFE (o) spectra by KOREL spectral disentangling, includ- 
ing LJ into the solution, versus orbital phase (Porb = 4^55970, 
last column Table 0}. Phase zero corresponds to the periastron 
passage. 
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Fig. 4. Normalised KOREL disentangled profiles of the Sim 
4553 A and 4568 A absorption lines of primary (black) and 
secondary (gray). 
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Fig. 5. The fflPPARCOS data folded according to the close or- 
bit of £ Lup (Porb = 4'?55970; zero phase corresponds to the 
epoch of periastron HJD 2447912). The maximum and mini- 
mum values are indicated by an asterisk. 
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Fig. 6. logTeff-log/^/i?© diagram for Schaller's et al. (1992, 
dashed fines) and Claret's (2004, fuU fine) evolutionary models 
with masses between 63Mq and 12.5M0, as indicated in the 
figure. The dotted line represents stars with an age of 15 Myr. 
The black box indicates the position of the primary of e Lup as- 
suming log Teff,i e [4.280; 4.364] K, and i?i e [3.6, 6.7]i?0. 
The dark gray box gives the restricted area assuming an age 
of approximately 15 Myr, while the figfit gray box is the cor- 
responding area of the secondary of e Lup. The two asterisks 
indicate the position of the representatives of the two compo- 
nents of £ Lup. 
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Fig. 7. Si III line-profiles observed at three different phases of 
the orbital period. The top (bottom) of each panel shows the 
'original' profile of the primary (secondary) in black and the 
'disentangled' primary (secondary) profile in gray. 
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Fig. 9. Lomb-Scargle periodograms of {v) of the disentangled 
Si III 4552.654 A profiles of the primary. The bottom panel 
gives the periodogram of the data after prewhitening with the 
highest peak at 0.06 c d~^. The top of each panel gives the 
window function, shifted towards the position of the highest 
frequency peak. 
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Fig. 10. Lomb-Scargle periodograms of the first normalised 
moment of the disentangled Si iii 4552.654 A profiles of the 
secondary. The moments were calculated using fixed integra- 
tion boundaries. The bottom panel gives the periodogram of 
the data after prewhitening with the highest peak at 0.28 c d~ ^ . 
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Fig. 11. The first normalised moment (dots) calculated from the 
Si III 4552.654 A profiles of the disentangled spectra of the sec- 
ondary, after prewhitening with peaks in the interval - 2 c 
d \ folded with 10.36 c d'^ . The full line is the fit with 10.36 
cd-i. 
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Table 6; RVs of the primary and secondary of e Lup. The columns give HJD, the 
RV of the primary, and an error estimate, the RV of the secondary and an error 
estimate, and a label indicating the source of the RV value. Error estimates of the 
RVs are only available for the recent data. As error estimate, we adopted the rms 
errors of the mean of the RVs calculated by CC from the Si ill, Mg ii, and He I 
profiles mentioned in the text. 



Date (HJD) 


RVl 


error 1 


RV2 


erTor2 


Note 


2417691.7678 


-37.6 




62.5 




A 


2417703.5991 


31.2 




-30.0 




A 


2418760.9114 


74.9 




-77.3 




A 


2418825.7282 


14.2 








A 


2418829.4971 


52.3 




-70.6 




A 


2419112.9000 


9.2 








A 


2419156.7000 


11.3 








A 


2419157.6720 


74.8 




-41.5 




A 


2419157.8070 


52.5 




-74.2 




A 


2419221.5900 


47.6 




-73.4 




A 


2435297.9640 


28.0 




-9.0 




B 


2436793.8760 


4.0 




-19.0 




B 


2436798.8560 


8.0 




-21.0 




B 


2438456.5951 


-42.0 




54.0 




C 


2438463.6087 


77.0 




-74.0 




C 


2438478.6309 


5.0 








C 


2438481.5881 


77.0 




-66.0 




C 


2438494.5070 


8.0 








C 


2438505.4036 


9.0 








C 


2438518.4670 


70.0 




-67.0 




C 


2438520.4831 


-28.0 




55.0 




C 


2438563.4078 


53.0 




-44.0 




C 


2438588.3454 


-12.0 




46.0 




C 


2438595.3089 


50.0 




-40.0 




C 


2438598.1947 


-33.0 




52.0 




C 


2438599.2336 


3.0 








C 


2438600.2796 


73.0 




-81.0 




C 


2438601.2175 


8.0 








C 


2438602.2524 


-21.0 




53.0 




C 


2438604.1792 


10.0 








C 


2438884.4871 


3.0 








C 


2438906.4533 


42.0 




-29.0 




C 


2439363.1908 


4.0 








C 


2439363.2438 


4.0 








C 


2439625.4602 


10.0 








C 


2439626.4213 


73.0 




-87.0 




C 


2439671.2942 


27.0 




-47.0 




c 


2439696.2405 


-3.0 








c 


2439696.3665 


-3.0 








c 


2439723.1773 


3.0 








c 


2439723.2742 


3.0 








c 


2439756.1914 


-33.0 




60.0 




c 


2440012.3611 


-31.0 




54.0 




c 


2440333.3966 


58.0 




-51.0 




c 


2440369.3663 


59.0 




-59.0 




c 


2449965.6037 


-38.4 


2.0 


47.7 


4.5 


D 


2450195.8947 


76.0 


2.1 


-66.1 


5.5 


D 


2450196.6035 


13.0 


3.5 


13.7 


6.5 


D 


2452770.5066 


-1.3 


8.2 






D 


2452770.5109 


-2.9 


7.7 






D 


2452770.5170 


-2.0 


7.6 






D 


2452770.5368 


-1.9 


7.3 






D 


2452770.5587 


-2.4 


6.8 






D 


2452770.6025 


-1.6 


6.2 






D 
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Table 6: continued 



Date (HJD) 


RVl 


error 1 


RV2 


error2 


Note 


2452770.6388 


-1.1 


4.4 






D 


2452770.6604 


-0.3 


4.2 






D 


2452770.6930 


-1.5 


4.2 






D 


2452770.7516 


-0.1 


4.6 






D 


2452770.7702 


-0.9 


2.2 






D 


2452770.7896 


-1.2 


2.3 






D 


2452770.8095 


0.0 


2.4 






D 


2452770.8296 


0.3 


2.6 






D 


2452770.8493 


0.7 


2.7 






D 


2452770.8704 


1.3 


2.7 






D 


2452770.8899 


0.5 


3.3 






D 


2452770.9119 


1.5 


2.8 






D 


2452771.5107 


46.7 


4.8 


-47.7 


5.7 


D 


2452771.5300 


48.7 


4.7 


-50.3 


7.5 


D 


2452771.5491 


50.5 


5.5 


-52.2 


8.6 


D 


2452771.5918 


53.3 


5.2 


-55.6 


11.3 


D 


2452771.6262 


55.4 


5.2 


-59.9 


17.4 


D 


2452771.6631 


59.0 


4.6 


-60.6 


7.8 


D 


2452771.7543 


64.6 


6.7 


-70.6 


11.6 


D 


2452771.7782 


66.3 


2.0 


-71.5 


12.4 


D 


2452771.8008 


66.8 


3.3 


-71.2 


7.1 


D 


2452771.8251 


67.8 


3.5 


-74.5 


7.8 


D 


2452771.8491 


69.2 


3.7 


-76.6 


11.7 


D 


2452771.8720 


70.6 


4.5 


-76.1 


13.0 


D 


2452771.8961 


70.5 


2.8 


-79.0 


13.3 


D 


2452772.6170 


11.7 


5.5 






D 


2452772.6505 


9.8 


3.4 






D 


2452772.7490 


5.0 


3.8 






D 


2452772.7740 


4.3 


3.7 






D 


2452772.8111 


2.1 


6.4 






D 


2452772.8455 


0.9 


5.4 






D 


2452773.6156 


-31.9 


4.8 


46.6 


9.1 


D 


2452773.6476 


-32.5 


4.8 


46.4 


9.1 


D 


2452773.6813 


-33.8 


4.6 


46.3 


9.6 


D 


2452773.7155 


-34.6 


3.2 


48.7 


9.6 


D 


2452773.7450 


-35.7 


3.1 


48.7 


9.3 


D 


2452773.7722 


-36.2 


3.2 


48.6 


9.1 


D 


2452773.7945 


-36.8 


3.3 


48.3 


7.4 


D 


2452773.8105 


-36.0 


4.3 


50.9 


7.2 


D 


2452773.8389 


-36.5 


4.5 


49.8 


7.4 


D 


2452773.8674 


-36.7 


4.7 


50.6 


11.3 


D 


2452774.5575 


-30.9 


4.1 


48.0 


14.7 


D 


2452774.5683 


-33.9 


2.1 


41.5 


8.4 


D 


2452775.5329 


3.6 


3.8 






D 


2452775.5555 


3.7 


3.7 






D 


2452775.5991 


5.2 


3.8 






D 


2452775.7024 


6.9 


2.6 






D 


2452775.7524 


8.3 


4.9 






D 


2452775.8086 


10.2 


13.4 






D 


2452775.8306 


11.5 


15.2 






D 


2452775.8556 


11.8 


15.8 






D 


2452780.2814 


1.5 


5.7 






E 


2452780.3087 


1.8 


14.9 






E 


2452780.3416 


2.2 


5.5 






E 


2452780.3824 


5.6 


8.8 






E 


2452780.5423 


32.8 


7.1 


-36.8 


14.4 


E 


2452780.5751 


34.8 


8.6 


-42.9 


5.7 


E 


2452781.3053 


61.3 


5.8 


-70.3 


11.0 


E 


2452781.3111 


60.4 


2.2 


-70.6 


5.9 


E 


2452781.3397 


56.3 


5.2 


-67.3 


13.3 


E 
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Table 6: continued 



Date (HJD) 


RVl 


error 1 


RV2 


error2 


Note 


2452781.3733 


54.7 


5.3 






E 


2452781.3930 


53.4 


6.0 






E 


2452781.4285 


50.0 


7.9 


-57.7 


8.9 


E 


2452781.4345 


51.1 


4.5 


-57.2 


8.2 


E 


2452781.4598 


48.8 


4.6 


-55.8 


7.3 


E 


2452781.4953 


43.2 


4.6 


-51.5 


4.9 


E 


2452781.5226 


40.4 


5.1 


-48.5 


3.0 


E 


2452781.5506 


38.5 


3.4 


-46.1 


4.1 


E 


2452781.5776 


33.8 


4.3 


-40.5 


8.9 


E 


2452781.6056 


32.6 


4.0 


-34.5 


12.1 


E 


2452783.2466 


-38.7 


4.8 






E 


2452783.2527 


-38.7 


4.6 






E 


2452783.2601 


-36.9 


3.7 


48.4 


13.8 


E 


2452783.2666 


-37.3 


3.6 


46.9 


10.6 


E 


2452783.2745 


-36.8 


2.7 


45.7 


8.9 


E 


2452783.2815 


-36.9 


2.6 


45.7 


9.3 


E 


2452783.3250 


-36.1 


3.2 


46.4 


3.9 


E 


2452783.3745 


-36.3 


2.9 


46.6 


11.2 


E 


2452783.4047 


-35.6 


3.2 


46.2 


4.8 


E 


2452783.4123 


-35.9 


3.0 


44.7 


6.5 


E 


2452783.4423 


-35.4 


1.3 


45.1 


7.3 


E 


2452783.4911 


-35.2 


1.8 


44.5 


7.1 


E 


2452783.4974 


-35.6 


2.4 


40.4 


5.9 


E 


2452783.5244 


-35.6 


1.8 


41.2 


6.2 


E 


2452783.5491 


-35.3 


1.2 


41.7 


7.0 


E 


2452783.5988 


-34.0 


6.0 


41.1 


13.9 


E 


2452784.2934 


-5.7 


5.4 


-6.8 


8.4 


E 


2452784.3026 


-5.3 


4.7 


-6.4 


8.2 


E 


2452784.4391 


-5.0 


5.4 


-2.1 


4.1 


E 


2452784.4463 


-4.8 


4.9 


-1.5 


3.3 


E 


2452784.4661 


-0.6 


4.2 


-2.6 


4.3 


E 


2452784.4708 


-2.2 


5.1 


-2.4 


4.5 


E 


Z-^-jZ / o'4-.'4- / o / 


-J.J 




-J.Z. 


J. 1 


Cj 


2452784.5072 


-4.5 


A.l 


-4.2 


5.0 


E 


2452784.5116 


-4.5 


4.9 


-3.0 


4.6 


E 


2452784.5481 


-4.3 


5.1 


-2.7 


5.5 


E 


2452784.5542 


-2.4 


3.5 


-2.6 


4.7 


E 


2452784.6059 


-1.9 


4.1 


-1.8 


5.1 


E 



A: CampbeU & Moore (1928); B: Buscombe & Moms (1960); C: Thackeray (1970); 
D: Average of the first moments of the Si iii 4553 A and 4568 A profiles; 
E: Median of the cross-correlation measurements calculated from the Si III 4568 A, 
Mg II 4481 A, He I 5016 A, He I 4917 A, He I 5876 A and He I 6678 A profiles 



